opar <- par(no.readonly = TRUE)
par(lty = 2, pch = 17)
drugA <- c(16, 20, 3, 45, 53)
plot(dose, drugA, type = "b")
par(opar)